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Abstract 

The electronic structure of nanowires in contact with metallic electrodes of experimentally rele- 
vant sizes is calculated by incorporating the electrostatic polarization potential into the atomistic 
single particle Schrodinger equation. We show that the presence of an electrode produces local- 
ized electron/hole states near the electrode, a phenomenon only exhibited in nanostructures and 
overlooked in the past. This phenomenon will have profound implications on electron transport 
in such nanosystems. We calculate several electrode/nanowire geometries, with varying contact 
depths and nanowire radii. We demonstrate the change in the band gap of up to 0.5 eV in 3 nm 
diameter CdSe nanowires and calculate the magnitude of the applied electric field necessary to 
overcome the localization. 

PACS numbers: 73.22.-f,73.40.Ns,73.43.Cd 



1 



Furture nanoelectronics will depend on electron/hole transport in a nanostructure, and 
crossing the nanostructure/metal electrode interface. Understanding the properties of such 
nanocontacts, especially new phenomena unique to the nanocontacts is thus of paramount 
importance to the future of nanoelectronics and it is currently an intensively studied sub- 
ject. Significant progress has been made in recent years in fabricating variety of nanostruc- 
tures and nanocontacts ranging from OD structures like isolated quantum dots, to more 
sophisticated structures, such as tetrapods, nanoribbons, nanorods and nanowires, etc. 

. Transport measurements in such nanostructures are often conducted by contacting 
a semiconductor nanostructure with large metallic electrodes, often deposited on top of the 



nanostructure 
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111 ] Hitherto, theoretical interpretation of such experi- 



ments has often been based on transport and electronic structure calculations with isolated 
nanostructures and ignoring influence of the electrodes. Even if the electrode is included in 
the calculation, its effects are often described by a short range correlation theory, such as 
local density approximation (LDA). For example, a nanocontact between an extremely small 
Si nanowire (with a diameter of 5 A) and Al electrode has been calculated by U. Landman 
et al. 12] . However, due to the use of LDA, the authors did not find the localization effects 
shown in this letter. The existence of the electrode will introduce the long range correlation 
effects well into the semiconductor nanowire, causing the change of its electronic structure. 
The long range correlation effects cannot be described by LDA, they can only be addressed 
by many-body theories, like the GW method, albeit the high cost of the GW method pre- 
vents it from being used to calculate the systems we discuss here. Fortunately, in the static 
approximation, the long range correlation effect allows a simple classical interpretation, [l^] 
It can be represented by an additional potential in the single particle Schrodinger equation, 
corresponding to a classical surface polarization potential (or image potential for metal). 
As we show below, this surface polarization potential is surprisingly large, particularly for 
small nanostructures, where it alters the familiar electronic structure and gives rise to a 
nanocontact phenomenon: electrode induced wavefunction localization. This localization 
disappears when the system becomes large, which explains why it has been overlooked in 
the past. 

In this letter we study the common case of the CdSe nanowire attached to a generic 
metallic contact. In practical calculations we use a long nanorod to represent a nanowire. 
The physical properties of CdSe nanorods such as optical spectroscopy, conductivity, electric 
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dipole, etc., have been extensively studied in the past decade Q, 
imentally and theoretically. The charging properties of isolated nanostructures (which are 
dependent on the surface polarization potential) have been studied theoretically, for quan- 



tum dots 
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21], and tetrapods [13J]. It has been demonstrated that the surface 



polarization potential (when the set up does not include a metallic electrode) plays an impor- 
tant part in the quantum dot charging energy. [l3] Here we study the electronic properties of 
the CdSe nanorods in contact with metallic electrodes taking into account surface polariza- 
tion potential. We use semi-empirical pseudopotential method (SEPM) to describe electron 
Schrodinger equation 22|. We assume that the metal electrode makes a Schottky contact 
with semiconductor nanorod as suggested in Ref. |23| (to avoid possible complications due 
to charge transfer). 

Within the SEPM we ordinarily solve the single-particle Schrodinger equation using 
planewave basis sets 

{- l fJ 2 + V}Ur) = EMr) (1) 

where, V is a potential which includes both local (sum of the screened atomic pseudopoten- 
tials) and nonlocal (s,p,d, and spin-orbit coupling) parts. The pseudopotentials are fitted 
to match the bulk experimental band structure. The potential V is bulk-like inside and 
zero outside the nanorod. The dangling bonds at the nanorod surface are passivated with 
a model ligand potentials to eliminate the midgap surface states. Once the potential V for 
the entire nanorod is constructed, the Eq. [1] can be used to solve for the eigenenergies and 
;he wavefunctions in the nanorod. In practice, we use the folded-spectrum method (FSM) 



24j in order to solve for only a few states in the valence and conduction band adjacent to 



the bandgap. The use of FSM allows us to perform calculations of large nanostructures 
containing tens of thousands of atoms, and approach experimentally relevant sizes range. 

In the presence of the electrode, the surface polarization potential due to interaction of 
an electron with its image charge becomes significant. Proper treatment of these long range 
correlations requires a many-body approach, such as GW method. However, it has been 
shown 13] that (under a static approximation) the self-energy potential in GW equation can 
be split into a bulk potential and a surface polarization potential. The GW equation is then 
split into two parts, the bulk Hamiltonian (analogous to EqJT]) and the surface polarization 
part, which is reduced to an electrostatic surface polarization potential P{r). Thus, the GW 
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equation can be approximated as 

{_I V 2 + V ± P(r)}ipi(r) = eMr) (2) 

where + and - are applied to electron and hole states, respectively, and are the quasiparticle 
energies. In the nanostructure the polarization potential P(r) can be derived from the many- 
body GW model, [3] which results in the same form as in the early electrostatic model of 
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P(r) = \ lim [W nano (r', r) - W bulk {r\ r)\ (3) 

where W nano (r', r) is the screened electrostatic potential in the nanorod at r' created by a test 
charge at r, and Wb u ik(r',r) is such potential in the bulk. In order to find the electrostatic 
potential W(r', r) we solve the Poisson equation for the structure shown in the inset to Figfj] 
(electrode size not to scale) 

V[e(r)V0(r)] = Anp(r) (4) 



. The die 
tructure 



ectric function e(r) is modeled following our study of dielectric response of a nanos- 
27], which equals ecaSe — 10 at the center of the nanorod, and approaches 1 near 



the surface of the rod. We use Dirichlet boundary conditions of 1/r at the domain bound- 
ary and and obtain P(r) which is then added to the previously computed CdSe nanorod 
potential V(r) in the Eq. [2} for the subsequent SEPM calculation of the electronic structure. 

We calculated surface polarization potential P(r) for nanorods length of 23 nm, and the 
diameters of 1.5, 3, 6, and 10 nm. The diameter of the electrode was much larger than 
that of the nanorod (<i e / ecir .=23 nm) in order to minimize the influence of a specific electrode 
geometry and mimic the situation of many experimental nanocontact setups. The SEPM 
calculations were performed for CdSe nanorods in the wurtzite crystal structure, length of 
23 nm, and 3 nm diameter, with total number of atoms of 5434. The polarization potential 
P(r) was generated for a) the nanorods embedded into a metal electrode by one half and one 
quarter of its length, b) nanorod and the electrode in contact, c) nanorod and the electrode 
separated by a 1 nm layer of vacuum, and d) a free standing nanorod. 

Figure [1] shows the polarization function P{ function of the ^-coordinate running; 

along the nanorod central axis, for nanorod diameter ranging from 1.5 to 10 nm, in the case 
of the nanorod embedded into the electrode by a quarter of its length. For the 3 nm rod we 
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FIG. 1: (color online) Polarization potential P(z) as a function of the z-coordinate running along 
the nanorod center axis (in case of the 3 nm rod, also the weighted average (P(z))). The nanorod 
is embedded into the electrode by a quarter of its length, nanorods of diameter 1.5, 3, 6, and 10 
nm are used. The inset shows an example of the modeled electrode and CdSe nanorod used in 
the present work (electrode size is not to scale, which has a diameter of 23 nm). The dielectric 
constant is e = 10 in the CdSe nanorod and e — * oo in the metal electrode. 

also computed the weighted average (P(z)) as 



where \ipi(x, y)\ 2 = J \ipi(x,y, z)\ 2 dz, and ipi(x,y,z) is taken either as a conduction band 
minimum (CBM) or valence band maximum (VBM) wavefunction (they do not make any 
practical difference). Thus, (P(z)) is a measure of the effective influence of the P(r) on the 
relevant wavefunctions. For the 3 nm rod the comparison of the central axis and weighted 
average polarization functions shows that the central axis P(z) is a good measure of the 
weighted average (P(r)). 

The large P(r) introduced by the presence of an electrode is evident quarter length along 
the rod. The electrode influence decreases with increased nanorod diameter. At 10 nm 
diameter a drop of P(r) near the electrode is almost indiscernible. This indicates that 
the localization effect which will be discussed below does not exist in a macroscopic bulk 
contact, which is why this has been overlooked before. However, for small nanorods the 
effect is surprisingly large, the P(r) provides a strong confining potential for electrons or 
holes (about 0.5 eV for 1.5 nm nanorod), and leads to an electrode induced electron/hole 
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FIG. 2: (color online) Contour plot of the lowest three CBM wavefunctions of the 3 nm diameter 
CdSe nanorod, (a) CBM wavefunctions in the free standing nanorod, (b) CBM wavefunctions in the 
nanorod with an electrode covering one quarter of the rod on the right side. The semitransparent 
rectangles indicate the position of the electrode. 

localization. 

Figure [2] shows the real space contour plots of the three wavefunctions adjacent to the 
bandgap in the 3 nm CdSe quantum rod conduction and valence band. In Fig. E(a) three 
wavefunctions are calculated for a free standing CdSe nanorod, while for Fig. Mjo) nanorod 
is embedded into the electrode (indicated by the semi-transparent rectangles) by one quarter 
of its length. In the latter case both CBM and VBM wavefunctions are localized by the 
polarization potential induced by the electrode. While for CBM states, the wavefunctions 
are just shrunk in the z-direction, for VBM states the second and third states show different 
nodal structure in comparison with their free rod counterparts, indicating strong state mixing 
and crossing. Due to the electron and hole state localization, the transport properties 
of small nanorods measured with attached electrodes will be different from familiar free 
electron-like picture, and more resemble the case of the Coulomb blockade. 

Another important consequence of the electrode induced polarization potential is the 
change in the value of the bandgap. Figure [3] shows the VBM and CBM states evolution 
as a function of the overlap between the nanorod and the electrode for the 3 nm diameter 
CdSe nanorod. In this case the value of the bandgap is reduced by approximately 0.5 
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FIG. 3: (color online) Evolution of the electronic states for the 3nm diameter CdSe nanorod as a 
function of the nanorod/electrode overlap, (a) CBM and VBM eigenvalues (bandgap), dashed lines 
indicate the bandgap computed ignoring surface polarization potential; (b) and (c) localization 
function (defined in the text) for the CBM and VBM states respectively. 

eV (FigJS^a)), from 2.78 eV to 2.26 eV as the rod is embedded into the electrode. The 
bangap value changes slightly when the nanorod and the electrode are separated by the 1 
nm layer of vacuum (E g = 2.74 eV). Once the electrode is in contact with the nanorod the 
change is more pronounced (E g = 2.64 eV) and becomes more significant with increasing 
the electrode/nanorod overlap, eventually saturating. In the absence of the electrode, the 
polarization potential P(r) of a free standing nanorod induces the bandgap increase of 0.54 
eV (from 2.24 eV to 2.78 eV) in comparison with the gap calculated ignoring the surface 
polarization (dashed lines in FigEJ). Thus, the existence of the electrode will also significantly 
affect the value of the quasiparticle band gap. Note that this quasiparticle bandgap equals 
the difference the electron affinity and ionization energy. It is different from the optical 
band gap. For optical band gap one has to include the electron-hole Coulomb interaction, 
which partially cancells out the polarization potential effects, especially for isolated spherical 
quantum dots. 

In order to quantify the electrode induced localization of the electron and hole wavefunc- 
tions we define a localization function 





where the wavefunction center of mass is z = J \tpi(x,y, z)\ 2 zdr. It shows how spread-out 
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FIG. 4: (color online) Weighted average of the polarization function (P(z)) for 3 nm diameter CdSe 
nanorod, half-covered by the electrode, plotted for several values of applied external electric field. 
The inset shows the real space contour of the CBM 3 wavefunction along with their eigenvalues 
for the electric fields of 0, 1.3, and 2.2 xlO V/m. The wavefunction is delocalized by the electric 
field, eigenvalue is shifted by 0.23 eV. 

the wavefunction is throughout the nanorod. This localization is presented in Figure [3] (b) 
and (c), for CBM and VBM respectively as a function of the electrode/nanorod overlap. 
It is interesting to note that localization of the wavefunctions is different from that of the 
free standing nanorod already when the nanorod and the electrode are 1 nm apart. The 
localization becomes significant for both CBM and VBM when the electrode is in contact 
with the nanorod. Throughout the electrode/nanorod overlap range the electron states are 
less localized than the hole states due to the differences in their effective masses. As the 
electrode covers more of the rod length, wavefunctions become less localized since the width 
of the P(z) confining potential is increasing. Thus there is an optimal electrode-nanorod 
overlap for maximally localized electron/hole wavefunctions. For the 3 nm nanorod this 
optimal overlap appears to be around 3 nm. 

In order to make the electron mobile again, one can apply an electric field along the wire 
axis to overcome the wavefunction localization. To estimate the magnitude of this electric 
field we apply a linear potential drop to the total potential V(r) in Eq. [2] between the ends 
of the nanorod, and calculate the resulting electronic properties. This is demonstrated in 
Figure H] where the weighted average of the polarization potential (P(z)) is plotted for 3 
nm diameter CdSe nanorod embedded into the electrode by half of its length. The field of 
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2.2 x 10 7 V/m (corresponding to a bias voltage of about 0.5 V applied across the 23 nm 
nanorod) is needed to delocalize the CBM wavefunction. The inset shows a representative 
(the third) CBM wavefunction contour plots for electric fields of 0, 1.3, and 2.2 xlO 7 V/m 
(bias voltages of 0, 0.3, and 0.5 V). The wavefunction in the latter case is dragged by the 
field towards the center of the rod, while its eigenvalue changes from -2.34 eV to -2.57 eV. 
Although it is no longer localized, it is still affected by the electrode, the structure of the 
wavefunction is not the same as in a free standing rod (compare to Fig. [2^). 

In conclusion, using atomistic pseudopotential method combined with electrostatic po- 
larization potential P(r) we have demonstrated the electrode induced localization of CBM 
and VBM states in CdSe nanorods (representing infinitely long nanowire). The effect is 
surprisingly large for small (1.5 - 3 nm diameter) nanorods but becomes insignificant as 
nanorod size grows (6 nm and up). The large polarization potential P(r) induced by the 
electrode also lead to the narrowing of the quasiparticle bandgap by 0.5 eV in 3 nm CdSe 
nanorods. We quantify the wavefunctions localization and compute the electric fields nec- 
essary to delocalize states near the CBM. We show that the mere presence of the electrode 
alters the nanorod electronic properties in comparison with a free standing nanorod, and 
this should be taken into account in interpreting the experimental transport measurements. 
This is a unique nanocontact phenomenon absent in its macroscopic counterpart. 

In order to experimentally confirm the effect of electrode induced localization we propose 
a measurement using one of the established wavefunction mapping techniques, such as scan- 
ning tunneling microscopy (STM), along with a conductivity measurements for a series of 
small nanorods (1.5 - 3 nm diameter), embedded in the electrods by different degrees. 

This work was supported by U.S. Department of Energy under Contract No.DE-AC02- 
05CH11231 and used the resources of the National Energy Research Scientific Computing 
Center (NERSC). 
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